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1.0  INTRODUCTION 

Wind  tunnel  tests  concerned  with  engine-airframe  integration  usually  use  a  cold  jet  to 
simulate  the  effect  the  real  hot  jet  has  on  the  airframe  pressure  distribution  and  resulting 
afterbody  pressure  drag.  Unfortunately,  there  is  no  reliable  method  available  for  predicting 
the  cold  jet  operating  condition  necessary  to  simulate  a  specific  hot  jet  operating  condition 
with  respect  to  afterbody  drag.  This  problem  has  been  studied  both  experimentally  and 
theoretically  for  many  years.  The  studies  have  been  primarily  concerned  with  axisymmetric 
configurations  as  shown  in  Fig.  1;  however,  it  is  expected  that  three-dimensional 
configurations  will  be  studied  in  the  future. 

The  experimental  results  presented  in  Refs.  1  and  2  clearly  show'  that  the  afterbody  drag 
obtained  with  a  hot  jet  is  significantly  less  than  that  obtained  with  a  cold  jet  operating  at  the 
same  pressure  ratio.  The  reason  for  this  difference  is  that  both  the  jet  plume  shape 
(blockage)  and  the  turbulent  mixing  (entrainment)  that  take  place  between  the  jet  flow  and 
the  external  stream  can  propagate  disturbances  upstream  through  the  subsonic  flow  region 
over  the  afterbody.  In  Refs.  1  and  2  it  is  shown  that  matching  the  initial  inclination  angle  of 
the  jet  plume  is  sufficient  to  correlate  the  hot  and  cold  jet  afterbody  drag  data  for  supersonic 
flow;  however,  for  subsonic  and  transonic  flow  an  additional  correction  for  the  entrainment 
effect  is  necessary.  In  Ref.  3  it  is  shown  that  cold  jet  afterbody  drag  for  an  annular  jet  and 
various  jet  nozzle  area  ratios  can  be  correlated  by  matching  the  maximum  jet  plume 
diameter  rather  than  the  initial  inclination  angle  of  the  plume.  However,  some  annular  jet 
configurations  appear  to  be  strongly  influenced  by  the  jet  entrainment  and  require  an 
additional  correction.  Experiments  showing  a  pure  entrainment  effect  on  afterbody  drag  are 
presented  in  Ref.  4.  This  was  accomplished  by  measuring  the  drag  obtained  with  both  a 
nitrogen  (N;)  jet  and  a  hydrogen  (H?)  jet  for  the  same  nozzle  and  afterbody  configuration. 
Because  N;  and  H;>  have  the  same  specific  heat  ratio,  the  plume  blockage  effect  is  identical 
for  equal  operating  pressures.  Thus,  the  only  effect  on  afterbody  drag  must  be  that  caused 
by  the  difference  in  jet  entrainment.  These  results  show'  that  decreasing  the  molecular  weight 
of  the  jet  causes  the  drag  to  decrease. 

The  specific  characteristic  of  the  jet-mixing  process  that  influences  afterbody  drag  is 
usually  thought  to  be  a  displacement  of  the  inviscid  jet  boundary,  which  changes  the  jet 
plume  blockage  (Ref.  1).  This  concept  was  utilized  in  the  viscous/inviscid  analytical 
technique,  presented  in  Ref.  5,  for  estimating  attached  nozzle-afterbody  flows.  For 
separated  nozzle-afterbody  flows,  the  Chapman-Korst  theory,  as  applied  in  Ref.  6,  suggests 
an  entrainment  parameter  based  on  the  pumping  capacity  of  the  jet  plume  which  involves 
not  only  the  rate  of  mixing  but  also  the  location  of  key  streamlines.  The  author  has  studied 
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this  approach  extensively  and  has  found  it  to  be  inadequate  for  subsonic  Mach  numbers.  In 
this  report  a  theoretical  estimate  is  presented  for  the  drag  caused  by  jet  entrainment,  based 
on  the  induced  velocity  produced  by  the  turbulent  mixing  process. 

2.0  ANALYSIS 


2.1  GENERAL  APPROACH 


It  is  well  known  (Ref.  7)  that  the  turbulent  mixing  of  a  jet  with  an  external  stream  will 
produce  a  small  induced  velocity  in  the  external  stream.  The  magnitude  of  the  induced 
velocity  has  been  estimated  in  Ref.  7  to  be  usually  less  than  5  percent  of  the  free-stream 
velocity.  Although  the  induced  velocity  is  small,  the  corresponding  induced  pressure  on  an 
afterbody  results  in  a  significant  change  in  the  afterbody  drag  coefficient.  A  rough 
theoretical  estimate  of  the  induced  drag  coefficient  can  be  obtained  on  the  basis  of  the 
following  assumptions: 


1.  The  induced  velocity  is  a  small  perturbation,  and  the  resulting  pressure  change 
can  be  computed  using  the  first-order  small  disturbance  equation. 

2.  The  external  flow  is  attached  to  both  the  afterbody  and  the  jet  plume,  and  the 
local  external  flow  velocity  is  equal  to  u^.  Variations  in  the  local  stream  velocity 
have  been  considered  and  found  to  have  a  negligible  effect. 

3.  Over  the  length  of  the  afterbody,  the  induced  velocity  is  assumed  constant. 

4.  The  induced  velocity  vector  is  perpendicular  to  the  average  jet  plume  boundary 
angle  between  the  jet  exit  and  maximum  jet  plume  position.  This  angle  is  taken 
to  be  one-half  the  initial  inviscid  jet  plume  angle. 

5.  The  afterbody  surface  is  conical,  as  shown  in  Fig.  1,  and  the  induced  pressure  on 
the  blunt  base  is  the  same  as  that  on  the  afterbody. 

6.  The  operating  conditions  of  the  jet  are  sufficient  to  produce  a  maximum  jet 
plume  diameter  that  is  equal  to  or  greater  than  the  jet  exit  diameter. 


The  induced  velocity  parallel  to  the  afterbody  surface  is,  by  geometry, 


(1) 
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The  minus  sign  results  from  taking  the  positive  direction  of  v  away  from  the  centerline  as 
shown  in  Fig.  1. 

v  Induced  Velocity 

\  Due  to  Mixing 


- ► 


Figure  1.  Typical  nozzle-afterbody  flow  field 

From  assumption  1  the  induced  pressure  coefficient  is 

An  -  6 


C  =  2 
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“ 

/  V  \ 

1 

/ 

(~ 

I  sin 

( 

f 
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Therefore,  the  induced  drag  coefficient  is 

r 


ACdb  =  -2 


1  - 


V, 


sin 


Ai/  +  $ 


N 


where 


*  =  — 


The  total  drag  coefficient  for  the  afterbody  can  be  written  as 


CDB1'  =  CDBO  T  ACDB 


(2) 


(3) 


(4) 


where  Cdbo  is  the  afterbody  drag  coefficent  for  <t>  =  0  and,  therefore,  is  only  a  function  of 
the  jet  plume  blockage.  It  should  be  noted  that  the  induced  drag  coefficient,  given  by  Eq. 
(3),  is  also  a  function  of  a  commonly  used  blockage  parameter,  i.e.,  the  initial  inclination 
angle  of  the  jet  plume,  (Av  +  0N). 
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2.2  INDUCED  VELOCITY 

Complex  flow  fields  involving  turbulent  mixing  such  as  the  afterbody  flow  field  shown  in 
Fig.  1  can  be  analyzed  by  superimposing  on  the  inviscid  flow  field  an  estimate  of  the 
turbulent  mixing  obtained  by  analyzing  a  much  simpler  flow  Held.  This  is  the  basis  of  the 
well-known  Chapman-Korst  component  analysis  as  applied  in  Refs.  5  and  6.  In  this  report 
the  simpler  flow  field  selected  for  the  turbulent  mixing  analysis,  shown  in  Fig.  2,  consists  of 
an  axisymmetric  jet  in  a  parallel-flowing  external  stream.  The  analytical  approach  is  the 
same  as  that  presented  in  Ref.  7  and  is  based  on  the  following  assumptions: 

1.  The  static  pressure  is  constant  along  the  control  volume  boundary. 

2.  The  initial  boundary  layers  are  neglected,  and  therefore  the  mixing  zone  velocity 
profiles  are  similar. 

3.  The  width  of  the  mixing  zone  is  small  compared  to  the  radius  of  the  jet. 

4.  The  Prandtl,  Lewis,  and  Schmidt  numbers  are  unity. 


Control  Vol.  I — ^  v  =nx)  /v 

|| — ►  Control  Vol.  II 

Free  Stream  jj  *  Uoo 

Ppo  =  Pn  (!—♦  _ 


v  =F1X)  Mixing  Induced  Velocity 


rL+ YLh 


rn  I1— 


UN  L  Mixing  Zone 


) _  trn+ 

r,f+  Yin 


Mixing  Velocity 
Profile 


L  Dividing  Streamline 

Figure  2.  Model  for  determining  the  induced  velocity  due 
to  mixing. 

The  conservation  equations  for  Control  Volume  1,  shown  in  Fig.  2,  are: 


Mass  Flow: 


PN  uN  An  +  =  ^77r«,  P»e  j*  vdX  -  Pn  uK  AU 

+  2n  S  Pur*  +  -  al) 

rll  v  ' 
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Momentum: 

PN  UR  “  Pc ou»(A~  -  An)  =  2nr<» f  vdX  +  PN  ul  Ajj 

•'o 

+  2n  J*  pu“tdr  +  ^  \oa  -  ALj  ^ 

V 

The  unknowns  in  Eqs.  (5)  and  (6)  are  ihe  induced  velocity,  v,  and  the  mixing  zone 
velocity  profile  shape  and  size.  The  mixing  velocity  profile  shape  was  determined,  as  in  Ref. 
7,  by  solving  a  highly  simplified  boundary-layer  equation. 

*  -  +  U  -  *J  t  0) 

where 

=  li  f  erf  T})  (8) 


and 


(9) 


Thus,  the  mixing  velocity  profile  is  expressed  in  terms  of  its  own  coordinate  system,  (X, Y), 
which  is  shifted  relative  to  the  lip  of  the  jet  by  the  distance  Ym,  Fig.  2.  The  nondimensiona! 
distance,  7jm,  is  determined  from  Eqs.  (5)  and  (6)  by  eliminating  the  induced  velocity  term. 
The  final  equation  is 


where 


(10) 
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Thus,  Tjm  is  independent  of  X,  as  are  all  the  properties  of  the  mixing  zone  when  expressed  in 
terms  of  rj.  The  actual  distance  Ym  is  determined  from  the  definition  of  7j  and  requires  a 
value  for  a,  the  mixing  similarity  parameter,  which  is  determined  in  Section  2.3.  In  Ref.  7  it 
is  shown  that  ?ju  =  3.0  and  rjL  =  -3.0  are  sufficient  for  r?m  to  be  independent  of  the  limits  of 
integration.  Evaluation  of  the  mass  flow  and  momentum  integrals,  (!])  and  (I2),  requires 
local  values  of  the  mixture  gas  constant,  Rg,  specific  heat  at  constant  pressure,  Cp,  and  the 
total  temperature,  Tt.  The  required  equations  were  derived  from  assumption  (4)  by  the 
method  presented  in  Ref.  6  and  are  as  follows: 


Jet  Species  Distribution: 


k  = 


Gas  Constant: 


Specific  Heat: 


Total  Temperature: 


R 


R 


—  =  k  +  (1  -  k) 

R  R 

SN  8n 


k  +  (1  -  k) 


c 

p 


oo 


c 

f>» 


vsK  / 


Substituting  Eq.  (5)  into  Eq.  (6)  yields 


01) 


(12) 


(13) 


(14) 


Differentiating  with  respect  to  X  yields 


-  C2 


(16) 
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Equation  (16)  shows  that  the  induced  velocity,  v,  is  independent  of  X  and  varies  inversely 
with  the  radius,  rM.  However,  in  this  report  only  the  induced  velocity  along  the  external 
stream  edge  of  the  mixing  zone  is  of  interest,  so  by  assumption  3,  r^  ~  rN,  and  Eq.  (16) 
becomes 


v 


u 


30 


7L 


+■ 


1 


(H) 


Because  the  initial  boundary  layers  have  been  neglected  in  deriving  Eq.  (17),  the  induced 
velocity,  predicted  by  Eq.  (17),  near  the  jet  exit  will  be  unrealistically  large.  To  correct  for 
this,  the  induced  velocity  predicted  by  Eq.  (17)  is  arbitrarily  reduced  by  a  factor  of  0.5.  The 
final  equation  for  the  induced  velocity  to  be  used  in  Eq.  (3)  for  estimating  the  induced  drag 
coefficient  is 


d> 


1 


2a 


(18) 


The  remaining  unknown  in  Eq.  (18)  is  the  similarity  parameter  for  turbulent  mixing,  a. 


2,3  TURBULENT  MIXING 


The  equation  for  estimating  the  similarity  parameter  for  turbulent  mixing,  a,  was  derived 
by  applying  the  method  of  Ref.  8  to  the  flow  field  shown  in  Fig.  2.  The  method  is  simply  an 
application  of  the  flow  conservation  laws  to  the  Control  Volume  II,  using  Prandtl’s  mixing 
length  hypothesis  to  represent  the  turbulent  shear  stress,  7-D,  acting  along  the  dividing 
streamline.  The  Conservation  equations  are  as  follows: 


Mass  Flow 

Pnc  ' 

Momentum: 


(A~  -  M 

=  2  rr  r  p 

co  r  rso 

C  v  dX  -  2tt- 

f  r  +  pM  UM  ( 

\  / 

4 

'  0  1 

%  X 

-  \ 


)  (19) 


Po c  -  xn)  +  2"  /  rDrndX  ^  277  f 


L 

77  g  p u  ^  rdr 

D 


H 


2*r»  u= 


/ 


ax 


(20) 
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The  location  of  the  diving  streamline  is  determined  by  the  following  equation  for 
conservation  of  mass  flow. 

(21) 

These  equations  are  solved  on  the  basis  of  the  four  assumptions  presented  in  Section  2.2. 
Applying  the  assumptions  to  Eq.  (21)  yields 


P N  UN  AN  =  P N 


u  N  AU 


2  n 


/ 


purdr 


(.,)  _  &H-: *■> 

V  %  V  Au  (i  _  c») 


It  is  important  to  note  that  tjd  is  independent  of  X;  therefore, 

R  T.  C 


/  v  n  1  u 

•  l  det. \  «d  ‘d  . 

V  dr>  L  R  T  C 

'  '  8  N  ciV  P  V 


are  also  independent  of  X. 


,x 


Eliminating  (  vdX  from  Eqs.  (19)  and  (20)  yeilds 


r&V'-eX' 


By  Prandtl’s  mixing  length  theory. 


Now  assume  that 


However,  by  definition, 
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and 


•■b-  7b»— 1 


'» ■  (K’0  (f) 


(22) 


(23) 


(24) 

(25) 
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(in 
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exa-Px 
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Equation  (24)  thus  becomes 

r» 

ph  4 


(27) 


(28) 


Because  the  fluid  dynamic  properties  along  the  dividing  streamline  are  independent  of  X,  rD, 
as  given  by  Eq.  (28),  is  also  independent  of  X.  Therefore, 


(29) 


Substituting  Eqs.  (28)  and  (29)  into  Eq.  (23)  yields 


The  remaining  unknown  in  Eq.  (30),  kijb,  is  determined  by  applying  to  Eq.  (30)  the  well- 
known  experimental  result  (Ref.  9)  that  a  =  12  for  planar,  single-stream,  two-dimensional, 
incompressible  mixing.  Therefore, 


3.0  COMPARISON  WITH  EXPERIMENT 


The  validity  of  the  theoretical  induced  drag  coefficient  cannot  be  evaluated  directly 
because  there  are  no  direct  measurements  of  this  drag  term.  As  a  result,  the  theoretical 
induced  drag  coefficient  was  combined  with  the  maximum  jet  plume  diameter  blockage 
parameter  of  Ref.  3  to  form  a  correlation  method  that  can  be  applied  to  a  wide  variety  of 
experimental  data  for  evaluation.  This  correlation  method  accounts  for  the  effects  of  jet 
area  ratio,  exit  angle,  total  temperature,  molecular  weight  and  ratio  of  specific  heats  on 


13 


AEDC-TR-79-85 


afterbody  drag  for  a  given  external  stream  Mach  number,  Reynolds  number,  and  afterbody- 
geometry.  In  this  report  correlation  of  the  effects  of  the  jet  thermodynamic  properties  is  the 
primary  concern. 

3.1  EFFECT  OF  MOLECULAR  WEIGHT 

The  correlation  method  was  applied  to  the  data  of  Ref.  4  to  predict  the  drag  of  an  H2  jet 
from  the  measured  drag  of  an  N2  jet.  The  predicted  results  are  compared  with  experiment  in 
Figs.  3a,  b,  and  c  for  the  15-deg  AGARD  afterbody  and  in  Figs.  4a,  b,  and  c  for  the  25-deg 
AGARD  afterbody  at  Mach  numbers  of  0.6,  0.9,  and  1.2,  respectively.  As  shown,  the 
correlation  method  w'orks  well  at  Mm  =  0.6  for  both  afterbody  configurations  but 
underpredicts  the  induced  drag  at  Mw  =  0.9  and  1.2  by  about  the  same  amount  for  each 
afterbody.  The  afterbody  pressure  distributions  presented  in  Ref.  4  indicate  that  the  flow  is 
attached  at  all  Mach  numbers  for  the  1 5-deg  afterbody  and  is  attached  only  aiM^  =  0.6  for 
the  25-deg  afterbody.  This  indicates  that  the  correlation  method  is  insensitive  to  the  type  of 
flow  over  the  afterbody,  which  is  a  desirable  characteristic  of  the  method. 

The  theory  indicates  that  the  drag  of  an  N2  jet  is  greater  than  that  of  an  H2  jet  because 
the  induced  velocity  for  the  N2  jet  is  toward  the  jet  plume  (negative)  and  for  the  H2  jet  is 
away  from  l he  jet  plume  (positive).  As  a  result,  the  induced  velocity  for  the  N2  jet  tends  to 
accelerate  the  flow  over  the  afterbody,  decreasing  the  pressure  and  increasing  the  drag, 
whereas  the  opposite  occurs  for  the  H2  jet.  It  is  important  to  note  that  the  direction  of  the 
induced  velocity,  as  given  by  Eq.  (18),  is  independent  of  the  similarity  parameter  for 
turbulent  mixing,  o,  since  a  is  always  positive.  Thus,  only  the  magnitude  of  the  induced 
velocity  is  determined  by  the  rate  of  mixing. 

3.2  EFFECT  OF  THE  RATIO  OF  SPECIFIC  HEATS 

The  correlation  method  was  applied  to  the  data  of  Ref.  4  lo  predict  the  drag  of  a  C2H4 
(ethylene)  jet  from  the  measured  drag  of  an  N2  jet.  The  predicted  results  are  compared  with 
experiment  for  both  afterbody  configurations  in  Figs.  3a  and  c  and  4a  and  c  for  M  =  0.6 
and  1.2,  respectively.  As  shown,  the  correlation  method  agrees  well  with  experiment  for 
both  afterbody  contigurations  and  Mach  numbers.  The  theoretical  induced  velocity  for  the 
C2H4  jet  was  found  to  be  essentially  equal  to  that  for  the  N2  jet;  thus  the  difference  in  drag 
between  the  two  jets  is  due  to  differences  in  jet  plume  blockage. 

Included  in  Fig.  3a  is  a  prediction  of  the  drag  of  an  H2  jet  using  the  measured  drag  of  a 
C2H4  jet.  As  shown,  this  predicted  drag  is  essentially  equal  lo  that  predicted  using  the 
measured  drag  of  an  N;  jet,  thus  verifying  that  the  correlation  method  is  applicable  for 
simuJtaneous  changes  in  both  molecular  weight  and  specific  heat  ratio. 
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3.3  EFFECT  OF  TOTAL  TEMPERATURE  AND  RATIO  OF  SPECIFIC  HEATS 

The  correlation  method  was  applied  to  the  data  of  Ref.  10  to  predict  the  drag  of  a  hot  jet 
from  the  measured  drag  of  a  cold  air  jet.  The  hot  jet  was  produced  by  burning  ethylene 
(C2H4)  with  air;  thus  both  the  temperature  and  ratio  of  specific  heats  were  different  from 
cold  air.  The  predicted  results  are  compared  with  experiment  as  solid  lines  in  Fig.  5  for  the 

C  •  EKpe1-! nnenf  ~  Ref  10 

- Esli-ated  Hoi  ‘rom  Cod  Data  »'th  Matcha:  t>p>DM 

- Estimated  Hot  Irom  Cold  Data  with  Matched  and  Entrainment 


b.  M  =  0.95 

o© 


c.  =  1.5 

Figure  5.  Comparison  with  experiment  of  estimated  drag  for  a  hot  jet 
from  cold  jet  data  for  the  15 -deg  AGARD  afterbody. 
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15-deg  AGARD  afterbody  and  in  Fig.  6  for  the  25-deg  AGARD  afterbody  at  Mach  numbers 
of  0.6,  0.95,  and  1.5,  respectively.  As  shown  for  both  afterbody  configurations,  the 
correlation  method  works  extremely  well  at  =  0.6,  underpredicts  the  induced  drag  at 
=  0.95,  and  overpredicts  the  induced  drag  at  =  1.5.  The  difference  between  the 
predicted  and  experimental  drag  is  partly  due  to  neglecting,  in  the  correlation  method,  the 
real  gas  effects  and  the  jet  total  temperature  gradient  which  occurred  in  the  experiment.  In 
addition,  the  correlation  method  does  not  account  for  spatial  variations  of  the  induced 


Enpenment'Sef  ID 


N  PR 

a.  M  =  0.6 


b.  =  0.95 


c.  M  =1.5 

oo 

Figure  6.  Comparison  with  experiment  of  estimated  drag  for  a  hot  jet 
from  cold  jet  data  for  the  25-deg  AGARD  afterbody. 
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velocity  which,  no  doubt,  become  important  at  Mo  =  1 .5.  In  spite  of  these  deficiencies,  the 
correlation  method  produces  a  significant  improvement  (except  at  M#  =  1.5)  over  merely 
using  a  correction  for  plume  blockage,  as  is  shown  by  the  dashed  line  in  Figs.  4  and  5. 

The  theory  predicts  that  the  induced  velocity  for  a  hot  jet  has  the  same  direction  as  that 
for  an  H2  jet,  which,  as  previously  discussed,  results  in  a  lower  drag.  This  suggests  that  a 
cold,  low  molecular  weight  gas  may  be  used  to  simulate  the  effects  of  a  hot,  relatively  high 
molecular  weight  gas.  A  brief  numerical  study  of  the  theoretical  relation  for  the  induced 
velocity,  Eq.  (18),  revealed  that,  for  a  given  7^,  the  induced  velocity  is  only  a  function  of  the 
product  RgN  T[N.  Since  the  induced  velocity  is  a  weak  function  of  7n,  then  on  the  basis  of  the 
assumptions  made,  the  two  major  parameters  that  affect  afterbody  body  drag  are  RgN  T,N 
and  Dp/D;m.  Of  course,  there  are  minor  effects  caused  by  9 N  and  Ap,  as  may  be  seen  from 
Eq.  (3). 

3.4  EFFECT  OF  AREA  RATIO,  EXIT  ANGLE,  AND  BLUNT  BASE  AREA 

The  purpose  of  this  application  of  the  correlation  method  is  to  illustrate  the  importance 
of  maintaining  the  external  afterbody  geometry.  The  correlation  method  w'as  applied  to  the 
data  of  Ref.  9  to  predict  the  drag  of  two  15-deg  AGARD  afterbody  configurations  (one  with 
and  one  without  a  blunt  base)  with  =  1.0  from  the  measured  drag  of  a  15-deg  AGARD 
afterbody  without  a  blunt  base  having  MN  =  1.83  and  ~  5  deg.  The  predicted  results  are 
compared  with  experiment  in  Fig.  7  for  Mw  =  0.6  and  1.5.  Only  one  theoretical  result  for 
both  afterbodies  is  shown  because  the  geometry  differences  are  not  sufficient  to  produce 
significantly  different  results.  The  predicted  drag  is  in  good  agreement  with  experiment  for 
the  afterbody  configuration  without  a  blunt  base,  which  is  the  same  configuration  used  to 
obtain  the  reference  drag  data.  The  blunt  base  significantly  increases  the  drag,  and  this 
cannot  be  accounted  for  by  the  present  correlation  method.  However,  the  correlation 
method  does  apply  to  configurations  having  a  blunt  base  since  all  of  the  previous  data 
presented  were  from  blunt  base  configurations.  Geometric  similarity  must  be  maintained  if 
the  correlation  method  is  to  be  valid. 

3.5  EVALUATION  OF  CORRELATION  METHOD 

The  improvement  realized  by  using  the  correlation  method  w'as  assessed  by  comparing 
the  results  with  the  two  commonly  used  experimental  correlation  methods  of  matching  NPR 
and  Dp/Dfg.  For  the  free-stream  Mach  numbers  considered,  the  accuracy  of  each  correlation 
method  was  determined  for  two  cases,  Dp/D^  »  1.0,  and  its  maximum  value.  The 
evaluations  were  made  at  Mach  numbers  of  0.6  and  1 .2  for  the  N;,  H2,  and  C2H4  data,  and 
at  0.6,  0.95,  and  1.5  for  the  hot  jet  data.  In  all,  15  evaluations  were  made  for  each  afterbody 
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configuration.  The  results  of  this  evaluation  are  presented  in  Table  1  for  each  afterbody 
configuration.  As  shown,  the  correlation  method  presented  in  this  report  (DP/DN  +  Entr)  is 
from  50  to  60  percent  more  accurate  than  matching  NPR  and  from  40  to  50  percent  more 
accurate  than  matching  Dp/DN;  also,  the  average  predicted  drag  is  within  10  percent  of  the 
correct  value. 

Experiment  -  Ref  10 
O  Mfj  ■  1  8325.  6N  ■  5  deg  (No  Basel 
2  •  1  0,  9fJ  ■  0  iiVith  Basel 

■  Mty  =1  0,  0N  ■  0  I  No  Basel 

- Esli mated  frljj  ■  1  0  with  Matched  Dp/D\j 


a.  M  =  0.6  b.  M  =  1.5 

90 


Figure  7.  Comparison  with  experiment  of  estimated  effects  on  drag  of  jet 
area  ratio  and  blunt  base  area  for  the  15-deg  AGARD  afterbody. 


Table  1.  Evaluation  of  Correlation  Methods 


I.  25° AGARD  Afterbody 

l.  Average  absolute  error,  AC^gr, 
over  all  conditions. 


?.  Percent  error  in  predlded 

II.  15° AGARD  Afterbody 

1.  Average  absolute  error,  ACDEfT, 
over  all  conditions. 


2.  Percent  error  In  predicted  CDBT 


Correlation  Method 

Average 

CDBT 

NPR 

DP,0N 

KrN 

0.0173 

0.0149 

0.0081 

0. 1287 

13.4% 

11.  ffb 

6.3% 

0.0134 

0.0120 

0.0056 

0.0603 

22.2% 

19.  % 

9.3% 
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4.0  SUMMARY  OF  RESULTS 

The  major  results  of  this  study  are  as  follows: 

1.  The  theoretical  induced  drag  coefficient  developed  in  this  report  provides  a  first- 
order  estimate  of  the  afterbody  drag  increments  due  to  differences  in  jet  mixing 
(entrainment)  caused  by  different  jet  properties. 

2.  The  correlation  method  used  is: 

a.  40  to  50  percent  more  accurate  than  the  correlation  method  based  entirely  on 
the  blockage  parameter, 

b.  accurate  to  within  10  percent  (average)  for  the  two  afterbody  configurations 
considered  in  the  Mach  number  range  from  0.6  to  1.5, 

c.  most  accurate  at  Mach  number  0.6  and  least  accurate  in  the  Mach  number 
range  from  0.9  to  1 ,2,  and 

d.  relatively  insensitive  to  the  type  of  flow  over  the  afterbody. 

3.  The  theory  indicates  that  the  drag  of  a  hot  jet  or  H2  jet  is  less  than  the  drag  of  a 
cold  air  or  N2  jet  because  the  jet  mixing  induced  velocity  is  away  from  the  jet 
plume  boundary  for  the  hot  and  H2  jets  and  is  toward  the  jet  plume  boundary 
for  the  cold  air  and  N2  jets.  As  a  result,  the  flow  over  the  afterbody  is 
decelerated  (decreased  drag)  by  the  jet  mixing  for  the  hot  and  H2  jets  and  is 
accelerated  (increased  drag)  by  the  jet  mixing  for  the  cold  air  and  N2  jets. 

4.  The  theoretical  induced  drag  is  a  strong  function  of  the  product  RgN  T,N  and  a 
weak  function  of  the  ratio  of  specific  heats, 
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NOMENCLATURE 

A  Area 

b  Width  of  mixing  zone 

C  Crocco  number 

ACdb  Induced  drag  coefficient,  i.e.,  entrainment  drag 

CDBo  Blockage  drag  coefficient,  i.e.,  0  =  0 
Cdbt  Total  afterbody  pressure  drag  coefficient  based  on  AB 
Cp  Specific  heat  at  constant  pressure 

D  Diameter 

K  Constant  in  Prandtl’s  mixing  length 

k  Mass  fraction 

(  Prandtl’s  mixing  length 

M  Mach  number 

NPR  Nozzle  pressure  ratio,  ptN/p ^ 

P,  Total  pressure 

p  Static  pressure 

Rg  Gas  constant 

r  Radius 

T,  Total  temperature 

u  Velocity 

u'  Small  perturbation  velocity  parallel  to  the  afterbody  surface 

v  Jet-mixing-induced  velocity  along  free-stream  edge  of  mixing  zone 

X  Distance  from  jet  exit  along  inviscid  jet  boundary 

Y  Mixing  zone  ordinate,  positive  toward  the  centerline 
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7  Ratio  of  specific  heats 

rj  Nondimensional  mizing  zone  ordinate,  cjY/X 

Nondimensional  mixing  zone  width,  ab/X 
jfni  i]  location  of  inviscid  jet  boundary 

0  Angle,  Fig.  1 

Ar  Prandtl-Meyer  expansion  angle  at  jet  exit 

q  Density 

a  Similarity  parameter  for  turbulent  mixing 

t  Shear  stress 

0  Nondimensional  velocity,  u/un 

4>  Nondimensional  induced  velocity,  Eq.  (18) 

SUBSCRIPTS 
B  Afterbody 

D  Dividing  streamline 

L  Low-speed  edge  of  mixing  zone 

N  Jet  or  nozzle 

P  Jet  plume 

U  High-speed  edge  of  mixing  zone 

oo  Free  stream 
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